Combinatorics of lattice animals by simple group— theoretic approach 
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Using the commutation relations in the locally free group we give a simple combinatoric derivation 
of the partition function of (l+l)-dimensional directed lattice animals with a fixed set of roots on a 
strip. Comparing our approach to the two-dimensional hard square model, we establish connections 
between the group-theoretic and the algebraic methods of the solution of some two-dimensional 
lattice spin systems. 
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I. INTRODUCTION 

The term "directed lattice animals" is used as a col- 
lective name for several related models describing the 
growth of aggregates, as the growth of molecular layers on 
substrates. Directed animals in the two-dimensional ge- 
ometry are structures of occupied and unoccupied nodes 
on a lattice strip of arbitrary length along one dimen- 
sion and with a finite width along the other one. It is 
convenient to align the lattice diagonally to the strip's 
edges. The occupied sites on the lowest row are the roots 
(the "source points") and the animal has to satisfy the 
condition that each occupied site can be reached from 
at least one root over a path which does only contain 
occupied sites and which never goes opposite to the pre- 
ferred direction. If we ascendingly label the nodes in the 
lowest row from left to right, the set of roots is given by 
the indices of the occupied nodes Q . Figure [T] shows an 
example. 




n = 3 

Figure 1: Directed animal on a strip of width n = 3 with set 
of roots C = {1, 2}. Occupied nodes are drawn black, the 
roots are highlighted grey. The vertical arrow indicates the 
preferred direction. 

In this work, we give a simple derivation of the parti- 
tion function of directed lattice animals with fixed root 
configuration on a strip of finite width. This problem 
has been solved exactly in Q by the diagonalization of a 
specially designed transfer matrix via algebraic methods 
borrowed from intcgrablc systems. Specifically, in 



the authors defined the transfer matrix T between row 
configurations represented in a spin space and computed 
the eigenvectors of T by finding the appropriate scat- 
tering i?-matrix as adopted in the theory of integrable 
systems. The result shows an exponential growth pro- 
portional to 3* in the limit of large numbers of occupied 
nodes s. 

The approach exploited in our work uses the corre- 
spondence between directed animals and heaps of pieces. 
These heaps have been introduced in Q as a geometric 
interpretation of the Cartier-Foata's partially commuta- 
tive monoid. The combinatorial solution described in our 
paper relies on the natural commutative structure of the 
"locally free" group |^], which describes available con- 
figurations of directed animals on a given lattice. The 
main idea consists of establishing a bijection between 
word counting in the locally free group and the particular 
configuration of the directed animal. The most involved 
part of the method is the correct account of the specific 
root configurations. 

By presenting our method we would like to direct at- 
tention towards possible hidden connections between the 
algebraic method developed in Q and our geometric 
group-theoretic approach based on the concept of the 
locally free group. 



II. THE MAIN RESULT 

We compute the partition function of directed animals 
with fixed roots by using the correspondence between 
lattice animals and a partially commutative semi-group, 
more precisely the locally free idempotent semi-group. 
We prove the following theorem: 

Theorem 1. The number Afn{C, s) of directed animals 
on a strip of width n with s occupied sites and set of roots 
C = {fci, . . . , fcf} with ki < k2 < • • • < ki is given by the 
expression 
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J\fJC,s) ^ } \ ' sin sin- TT sin ^ ' ^ — 1 + 2 cos (1 
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where Hq ~ ki ~ 1, /i^ = k^^i — hi — 1 for 1 < i < I and 
hg = n — ke + ^. 

This paper is organized as follows. In section lllli we 
recall the correspondence between directed animals and 
heaps of pieces and give the definition of the locally free 
group. In section lTvl we consider lattice animals with a 
single root before generalizing to arbitrary numbers of 
roots in sectionlVl Finally, we summarize our results in 
section lVTl 

III. GROUP-THEORETIC FORMULATION OF 
THE PROBLEM 

A lattice animal as shown in figure [1] can equivalently 
be drawn as a heap of 2-dimensional rectangular boxes 
which have discrete horizontal positions (columns) and 
which can only touch each other at their corners. The 
roots of the heap are the boxes which stand on the bot- 
tom of a column. Figure [2] shows the equivalence of the 
two pictures. 

We write a;„({C}, s) for the number of heaps of s boxes 
on n columns with a set of roots C. As can be deduced 
from figure [2l its correspondence to the partition func- 
tions of directed animals is 



AA„({fci, . . . }, S) = UJ2ni{2ki - 1, . . . }, s) (2) 




1 2 3 4 5 6 



Figure 2: A lattice animal on a strip of width 3 drawn as a 
heap of boxes on 6 columns. 

The combinatorial treatment of heaps of pieces was 
formalized by G. Viennot The boxes which we con- 
sider (called dimers in are a special kind of pieces. 
Viennot showed that each heap of pieces can be seen as 
the graphical representation of an element of a partially 
commutative semi-group. For our heaps of boxes, this is 
the locally free idempotent semi-group. 



I 

Definition 1 (Locally free idempotent semi-group). 

A semi-group of rank n with the generating set 
{/i, . . . , /„} is called a locally free idempotent semi-group 
CJ-I^_^_i if it satisfies following relations. 

f.fj = for >2 (3a) 

{f^? = for i^l...n (3b) 

and if each two generators fi, fi+i, i = I . . . n—1 produce 
a free sub-semi-group of CTI^^t^. 

We will always refer to the generating set specified in the 
definition above. Each of the generators corresponds to 
the act of dropping a box in one of the columns. Figure [3] 
shows the meaning of the commutation relations in the 
heaps picture. 




Figure 3: The uppermost picture shows that the dropping 
of boxes on adjacent columns does not commute, whereas it 
commutes if there is at least one column inbetween (middle) . 
The sequential dropping of two boxes in the same column is 
the same as dropping only one box (below). 

Each element of CTT^j^^ can be written as a product 
of the generators like 

9 ~ fai ' /q!2 ' ' ' fas (4) 

Here, s can be any natural number. Each such product 
can be written as a word 

W = {ai,a2,...,as) (5) 

However, due to the partial commutativity of the gen- 
erators, the word representation of a semi-group element 
is in general not unique. Two words which represent the 
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same semi-group element belong to the same equivalence 
class. For example, it is 



(1,3,2)^(3,1,2) 



(6) 



In order to find a unique word representation of each 
element of £J^I^_|_i, the authors of ^ introduced the or- 
dered word form. 

Definition 2. An s~tuple W = {ai^ . . . ,as) with ai E 
{I, . . . ,n} for i = I . . . s has the ordered form, if and only 
if it satisfies the following conditions. 



if ai — 1 then a^+i G {2, . . . , n} (i) 

if 1 < ai < n then aij^i G 

{ofj - l,ai + 1, . . . ,n} (ii) 
if ai = n then a^+i = n — 1 (iii) 



e„(x,A) = (A + 1) — 



sm 0n 



sin (/){n + 2) 



In particular, 



e„(n. A) 



sm 0n 



sin (/){n + 2) 



(9) 



(10) 



Here and in the following, the angle cj) is defined via the 
identity 



arccos ■ 



VA+T 



(11) 



In eqs. ([5|- pII)) A G [^1,3] and all poles lie within this 
interval (8|. 

In order to derive results for the lattice animal problem 
it is necessary to find out how a given root translates 
to the corresponding word representation. This rather 
nontrivial task will be adressed in the next two sections. 



The conditions (i)-(iii) are visualized in figured 




1 ... x-1 X x+1 ... n 



Figure 4: The arrows point to the values the (i + l)th letter 
of the ordered word can take if the ith letter is x. 

Each element of CJ-X^^i is represented in one and only 
one way as an ordered word. The formal expression for 
the number 9n{x,s) of words of length s with first letter 
X in the box of size n is 

0n{x,s) = #{{x,a2,...,as) I (i)-(iii)} (7) 

The function 0„(x, s) satisfies the recursion relation 
which straightforwardly follows from the conditions (i)- 
(iii) and figure S) 

9n{x,S +l) = 0n{x- 1,S) +J2y=x+l^n{y,s), 
^„(a:,l) =4,.o: (8) 

6»„(0,s) =e{n + l,s) =0. 
The generating function 

oo 

e„(x,A) =^^^„(a;,s)A-^ 

s=l 

can be easily computed by solving ^ (see [1] for details) . 
It reads 



IV. HEAPS WITH A SINGLE ROOT 

In this chapter we search for the number of heaps with 
a single root at a given site. The question we ask can 
be rephrased as follows: "Which condition (additional 
to conditions (i)-(iii)) does the word representing a heap 
with one single root satisfy?" 

Suppose we want to generate a heap with a single root 
in the fcth column. Obviously, the first letter of the word 
has to be k. Since we do not want that the second box 
falls on the ground, we have to make sure that it falls 
ontop of the first one. Consequently, 012 has to be either 
k ^ I or k + I. In order to prevent that 012 takes any 
other value than one of these two, we put a "wall" at site 
k + 2 (see figure [5]) . As soon as the first box is dropped 
in column k + 1, we can shift the wall to the right by 
one step, because then, any box dropped in the k + 2th 
column will fall ontop of the one in column fc + 1 . When 
the first box is dropped in column k + 2, we can shift our 
wall again, and so on. 




1 ... x-1 X x+1 x+2 ... n 



Figure 5: Reduced set of possible values of Q2 compatible 
with the ordered word form and preventing of having roots 
other than in position qi = x. Jumps from position ai = x 
to positions a2 — x + 2, ...,n are blocked. 

Algebraically, the condition which ensures that the or- 
dered word represents a heap with a single root located 
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in the column x (see figure [5]) , reads 

aj < max {ai\l < j < i} + 2 V j G {2, . . . , s} (iv) 



Condition (iv) is automatically satisfied for ai = n. Con- 
sequently it is 



uj„{{n},s) = 6'„(n, s) 



(13) 



With (iv), we can now write the partition function of 
heaps in n columns with s elements and a single root in With the use of equation dTS]), the following theorem can 
the xth column (see figure [S]) as be derived. 

Theorem 2. The partition function ujn{{k} , s) of heaps 
ujn{{x}, s) — ^ { {x, . . . , Us) I (i)-(iv) } (12) with a single root at site k satisfies the equation 



w„({fc}, s) = 0fc(fc, s) + ^ 0k{k,so)9k+i{k + l,si 



(14) 



so+si=s 



soH hs„-fc=s 



r 



Since the proof is straightforward, it is not carried out 
here. 

The series (fT4|) has the form of a convolution, so we 
step over to the generating function 



r!„({fc},A) = ^c^„({fc},s)A-^ 

s=l 

for which equation (1141) transforms to 

n — k d 
d=0 j=0 

By inserting we get after some algebra 
sin k(f> sin(n — k + 



(15) 



(16) 



nni{k},\) = 



sin0 sin(n + 2)0 



(17) 



for A e [-1,3] 



V. GENERALIZATION TO ARBITRARY 
NUMBERS OF ROOTS 

Now we are in the position to generalize our results to 
arbitrary numbers of roots. 

Suppose that a heap of s pieces possesses two roots in 
ki and fc2 and that ki < k2- The representing word has 
the form 

W = (fci, a2, . . . ,asi, ^2, asi+2, ■ • • ,^31+82) (18) 

with si + S2 — s and 7^ fc2 V i = 1 . . . si. Since the 
first part ofW, (/ci , a2 , . . . , ) , does not generate a root 
except in fci, it has to satisfy condition (iv). On the other 
hand, agj+i = k2 generates a root, therefore it has to be 

max {ctj \j < si} < k2 — 1 



Similarly, the second part, (A:2, must not gener- 

ate another root except in k2 and therefore it has to sat- 
isfy (iv) as well. Consequently, the partition function for 
two roots has the form 

U!n{{kl,k2}, s) = 

^ '^/c2-2({fcl},Sl) W„({/C2},S2) (19) 

Sl+S2 = S 

The generating function reads 

a.({fci, fc2}, A) = r!fc,_2({fci}, A) r!„({fc2}, A) (20) 

In this way, we can proceed to heaps with arbi- 
trarily many roots. The generating function of heaps 
with £ roots ki < k2 < ■ ■ ■ < kg reads 

VLn{{ki,k2,...,ki},\) = 

f^n({fca,A)n"fc.+i-2({fcJ,A) (21) 

i=l 

Inserting pT|) . and setting 

Ao = fci (22a) 
Ai = h+i -h-l for 1 <i <l (22b) 
Ai = n-ke + 1 (22c) 



we obtain 
^n{{ki, 



,ki},X) = - 



sin(n + 2)(j) -^-^ sm 



nsinA,-^' 
-— F (23) 



for A £ [—1,3]. For £ = 1, equation ([23| reproduces 
equation ([TT]). 

The final step is to apply the inverse Z-transform to 
the generating function (j23l) in order to obtain the par- 
tition function. This is done by performing the integral 



Wn(C,s) = — 



nnic, A) 



27rz./];,|=p A--+1 



dX 



(24) 
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where the contour |A| = p surrounds all poles 



Am = 4 cos 



n + 2 



1 ; m = 1 



n+1 



The residues can be easily evaluated, and the partition 
function of heaps on n columns with set of roots C = 
{ki, . . . ,ki} is derived as 



(25) 



of f7„ (C, A) . Applying the residue theorem, we obtain 



=^Res 



sin A, 



sin(n + 2)0 sine 



n 



(26) 



--^^^^^=T.hT^''''^2h^2) n l + 2cos-^ (27) 



m—l 



Combining (|27|) and ([2]) , we arrive at our final result ((TJ . 
It coincides with the one given in [2j. 



VI. SUMMARY 

The number of directed animals with s occupied sites 
on a strip of width n for given root C has been calculated 
by using the analogy to the model of heaps of pieces. The 
result, which is formulated in theorem [TJ coincides with 
the one obtained by Hakim and Nadal [2] who used an 
algebraic method dealing with the diagonalization of the 
transfer matrix for some spin system. 

It has been shown by D. Dhar in that the statis- 
tical weight of a configuration of two-dimensional di- 
rected animals is the same as the statistical weight of 
a spin configuration of an antiferromagnetic Ising model 
with anisotropic nearest-neighbor coupling in an external 
magnetic field on a triangular lattice with a Hamiltonian 



(28) 



The correspondence between directed animals and the 
2D Ising spin system is set as follows. Define the fugac- 
ity, z, conjugated to the animal's size n (z is the same as 
A""'^ in our computations). Now the generating function, 
A{z) of directed-site lattice animals on the square lattice 
reads A{z) — —daf{a, 2:)|q=o where ln(l + z) + f{a, z) is 
the Gibbs free energy per site for the antiferromagnetic 
Ising system with the Hamiltonian H defined above. The 
parameters a and z are connected to the coupling con- 
stant L and to the chemical potential /i in (|28p by the 



equations e — 1 -I- z. 



-(l-a)z(l + z)- 



Note that A{z) coincides up to the redefinition z — >■ 
A-i with e„(n,A) defined in (IlOl)-(ITl|). The minor dif- 
ference between A{z) and 8„(n, A) consists in the fact 
that A{z) is defined for periodic boundary conditions, 
while 6„(n, A) - for open ones. 



By using the concept of the locally free group, the enu- 
meration of directed lattice animals becomes a rather 
simple combinatoric problem. So, it seems an intrigu- 
ing question how far the ideas of geometric group theory 
could be extended in the area of two-dimensional exactly 
solvable spin systems. 

The proposed combinatorial approach to the statis- 
tics of lattice animals admits few possible generalizations. 
The simplest one consists in relaxing the idempotent con- 
dition ff = fi (see I3bp , or replacing it by more general 
ones: either ff = or ff — 1 {p — 1, 2, ...). Such a modi- 
fication of the problem allows to consider directed lattice 
animals on a class of decorated lattices. For example, 
skipping the relation ff — fi, we pass to the triangular 
lattice - see Fig. [B^. 




Wn 



(b) 



)^ ^ )^ 



Figure 6: (a) Animals on triangular lattice; (b) Left-, Right- 
and Double-handed supports of particles on a square lattice. 

A more involved generalization of the model allows 
to distinguish between "left-handed" (L) and "right- 
handed" (R) supports of deposited particles as shown 
in Fig. To proceed with the description, suppose 

that at the initial time t = the system is empty. Then, 
at each time step, t — 1,2, . . . , N , an elementary unit 
("particle") is deposited at some column. Consider the 
deposited particles as distinguishable (i.e. "colored"), 
and assign to them the variable t according to the natu- 
ral time ordering. Now each particle carries two indices: 
the index t fixes its occurrence in time (i.e. a "color"), 
and the index ai labels the column in which the parti- 
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cle is deposited. Now proceed with normal ordering, but 
apply it twice to the sequence: first - exactly as it is de- 
scribed in Definition 2, i.e. pushing the smaller indices 
as left as possible (allowed by the commutation relations) 
and second - inversely, pushing the smaller indices ai as 
right as possible (also consistent with the commutation 
relations). The first ordering will catch the L-supported 
particles, while the second ordering - the R-supported 
particles. Since the particles are colored, we can easily 



identify whether a given particle has L, R, or double sup- 
ports. Introducing different weights WL,ifR,WD for L-, 
R- and double supports, we could distinguish between 
these configurations, which is crucial for a wide class of 
models (for example, for cellular automata ^6]). 

The authors are grateful to A. Vershik and R. Bikbov 
for helpful discussions. N.H. acknowledges the LPTMS, 
S.N is partially supported by the ANR grant 2011-BS04- 
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